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Abstract 

The measurement of the (3 asymmetry parameter in nuclear (3 decay is a potentially 
very sensitive tool to search for non V — A components in the charge-changing weak 
interaction. To reach the required precision (percent level) all effects that modify 
the emission pattern of the (3 radiation, i.e. the geometry of the setup, the effect of 
the magnetic field on the trajectories of j3 particles as well as (back) scattering in 
the source, on the sample holder and on the detector, have to be correctly taken 
into account in the analysis of the data. A thorough study of these effects and a new 
method based on detailed GEANT4 Monte-Carlo simulations that was developed for 
this purpose is presented here. The code was developed for j3 asymmetry measure- 
ments by means of the Low Temperature Nuclear Orientation (LTNO) method, but 
can in principle be generalized to other experimental setups using other polarization 
techniques. 
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1 Introduction and motivation 



In nuclear ft decay the search for physics beyond the Standard Electroweak 
Model is a key issue in a series ofprecision experiments that are carried out 



or being prepared now (see e.g. One of the powerful techniques to 



probe the structure of the Standard Model is a measurement of the angular 
distribution of ft radiation with the Low Temperature Nuclear Orientation 
(LTNO) method. One of the best examples of this is the famous experiment 



carried out by C.S. Wu et al. in 1957 (1121 ) and later repeated with better pre 



cision by Chirovsky et al. (Il3l ). The observed angular distribution was found 
to be in agreement with the expected 1 + AP (v / c) cos (8) dependence from 
the V-A theory of the weak interaction (the explanation of the parameters 
in this formula is given below). Such a measurement of the angular distribu- 
tion is a potentially very sensitive tool to search for non V — A (i.e. tensor 
and scalar) components in the weak interaction. However, present upper lim- 
its for such contributions (which are at the level of 8% (95 % C.L.) in the 
amplitudes (14)) require the angular distribution coefficient A (see below) to 
be determined with a precision better than about 1 %. In the past we have 
performed already several measurements of the angular distribution of ft radi- 



ation with the LTNO method (1171 ). (I la ). (Il9l ). However, in order to reach the 



1 % precision level the method had to be further optimized. Therefore, a new 
LTNO setup was recently installed at the Katholieke Universiteit Leuven (flil ). 
while a GEANT4 based Monte Carlo code for this type of measurements was 
developed as well. These Monte-Carlo simulations turned out to be a crucial 
factor, providing a detailed understanding of systematic effects such as e.g. 
scattering and magnetic field related effects. The method presented here is 
primarily meant to search for a tensor current admixture in the V-A structure 
of the weak interaction. It is also being used for the analysis and interpreta- 
tion of measurements performed at the ISOLDE-facility at CERN (jig) and 
can in principle be generalized to other experimental setups using different 
polarization techniques as well. 

The measured angular distribution of ft radiation emitted by oriented nuclei 
can be expressed as (jl5l ): 



W(6)= ^ )po ' =l + V -A.p.Q 1 .cose (1) 

N (0) 'unpol C 

Here N {0) pol (N {0) unpo] ) is the count rate in the detector at an angle 6 with 
respect to the nuclear orientation axis for polarized (unpolarized) nuclei, v/cis 
the ft particle velocity relative to the speed of light, P is the degree of nuclear 
polarization and A is the angular distribution coefficient that is sensitive to 
the possible presence of tensor currents. The modification of the angular dis- 
tribution pattern of the ft radiation due to the finite sizes of the detector and 
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the source, (back) scattering of the (3 particles and the presence of a magnetic 
field is taken into account by the so-called solid angle correction factor Qi. 
However, as will be pointed out below, the general approach to calculate this 
Qi factor (1201 ) is not applicable anymore when the influence of the disturbing 
effects becomes significant and/or high precision is required. Therefore, a new 
method based on Monte-Carlo simulations with the GEANT4 toolkit (12 ll ) was 
developed to correct the anglular distribution for the above mentioned effects. 
As the first measurements of the j3 asymmetry parameter with our new LTNO 
setup were performed with 60 Co, most of the simulations were done for this 
isotope. The experimental setup was described previously already in detail 
(flol). Only a short description will therefore be given here. This paper will 
deal mainly with the GEANT4 simulations. The performance of GEANT4 
for electron backscattering and electron energy deposition in the energ y ra nge 
from 0.040 MeV to 1.0 MeV was recently investigated by Martin et al. Q O) 
and by Kadri et al. (24). Good agreement with experimental data at the level 
of a few percent was obtained in all cases. Recently, GEANT4 simulations 
were also performed for scintillation detectors, leading to similar conclusions 



2 Experimental setup and solid angle correction factor 



In order to orient nuclei the Brute Force LTNO method is used here. The setup 
consists of a 3 He- 4 He dilution refrigerator, which permits to cool samples to a 
temperature of about 5-10 mK, in combination with a superconducting mag- 
net. A schematic view of the lower part of the setup is presented in Fig. [TJ 
The magnet was manufactured by Oxford Instruments Ltd. and can produce 
a field up to 17 T. (3 particles emitted by the radioactive nuclei are focused by 
the strong magnetic field and registered by the particle detector. The choice 
of particle detector is determined by its ability to work both at temperatures 
close to liquid He temperature and in a strong magnetic field. Si PIN photo- 
diodes produced by Hamamatsu Photonics (type S3590-06) with a thickness 
of 500 fim and an active surface of 9x9 mm 2 have been tested and showed 
good behavior under such conditions (}26l ). Detectors of this thickness are able 
to fully stop electrons with energies up to about 350 keV. Note that the 500 
/im thickness is thick enough for the isotopes of interest for this type of mea- 
surements (e.g 60 Co, 95 Nb, 133 Xe, etc.) and at the same time thin enough to 
have low sensitivity to 7 quanta. The source-to-detector distance is about 20 
cm and with the central field of 17 T the magnetic field at the site of the Si 
detector is about 0.8 T. 

The sample for measurements of the 60 Co f3 asymmetry parameter was pre- 
pared by diffusing 60 Co activity into a copper foil with a thickness of 20 ^m 
and a purity of 99.99+ %. A special temperature regime was used to limit 
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Fig. 1. Schematic view of the lower part of the experimental setup, showing also the 
detailed geometry of the sample holder and the position of the source foil (containing 
the 60 Co and 57 Co activities) and of the 54 Mn thermometer. The magnetic field is 
pointing vertically. 



the depth of diffusion to a few microns (1271 ) . This foil was then soldered onto 
the oxygen-free copper sample holder which is, via the cold finger, in thermal 
contact with the 3 He- 4 He mixing chamber where the lowest temperature of 
the refrigerator is reached. Geometrically the foil is positioned in the center 
of the magnetic field and perpendicular to the field direction (Fig. [T]). 

To determine temperatures in the millikeivin region a nuclear orientation ther- 
mometer (e.g. 57 Co(Fe) or 54 Mn(Ni)) is soldered on to the sample holder. The 
temperature is then deduced from the degree of anisotropy observed for the 
7 rays of this thermometer. The 7 rays are registered by a high-purity Ge 7 
detector installed outside the refrigerator at a distance of about 30 cm from 
the center of the magnetic field (see Fig. [I]). 

From the above given details of the experimental setup it can be seen that the 
following effects can modify the angular distribution pattern of the (3 radiation: 

- scattering of the f3 particles in the host foil material; 

- backscattering on the layer of solder; 

- backscattering on the sample holder; 
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- backscattering on the detector; 

- magnetic field influence on the (3 particle propagation; 

- possible Compton background from 7 rays accompanying the (3 decay and/or 
7 rays emitted by the nuclear orientation thermometer. 

As was already mentioned, the distortion of the emitted angular distribution 
of the (3 radiation is taken into account by introducing the solid angle cor- 
rection factor Qi (see Eq. [1]). A general formalism to calculate the Q\ factor 



is described in Ref. (I20l). For an axially symmetric detection system this Q± 
factor is defined as 

_ /sin(fl)cos(fl) £ (fl)dfl 
Ql ~ J sin (9) e{9) [Z) 
where 9 is the angle of the j3 particle emission direction with respect to the 
orientation axis (which coincides in this case with the detector axis) and e (9) 
is the detector response function which describes the efficiency of the detector 
for radiation incident at angle 9. The problem of determining Q\ thus reduces 
essentially to the knowledge of the detector response function for the j3 radia- 
tion. However in the presence of (back) scattering and a strong magnetic field 
this approach is not valid for precision measurements. A different and more 
inclusive method for correcting the angular distribution function was therefore 
developed. 

The new method consists of a Monte-Carlo simulation of the complete exper- 
iment. The output of this provides the (3 spectrum registered in the (3 particle 
detector for both isotropic emission of the (3 radiation (i.e. no nuclear po- 
larization) and anisotropic emission (i.e. when the sample is polarized). The 
angular distribution of the emitted (3 radiation is calculated for exactly the 
same degree of nuclear polarization P (see Eq. [I]) as obtained in the exper- 
iment (which can be deduced from the observed anisotropics of the 60 Co 7 
rays) and assuming the Standard Model value for the (3 asymmetry parame- 
ter, A SM . The simulated isotropic and anisotropic spectra are then analyzed 
in exactly the same way as the experimental spectra and the (3 asymmetry 
parameter A is extracted using the following relation: 



A exp l-W{9) 



exp 



Asm l-W{9) 



sim ' 



(3) 



where W (9) exp and W (9) sim are the anisotropy functions (Eq. [T]) obtained 
experimentally and from the simulations, respectively. 



To perform these simulations the GEANT4 toolkit (12 ll ). which was designed 
for the simulation of the passage of particles through matter, was used. The 
development of the simulation code was performed in several steps. Firstly, 
a simple simulation routine to calculate electron backscattering coefficients 
was developed. With this routine the influence of the different GEANT4 pa- 
rameters, i.e. the cut-for-secondaries (CFS) parameter and the so-called f r 
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parameter, on the calculated backscattering fractions was studied. The CFS 
parameter determines the production threshold for secondary particles, while 
the f r parameter limits the step size for tracking of electrons at the bound- 
ary between two materials. These parameters were then fixed at the values 
that provide the best match with the reference data that were obtained from 
the literature. Secondly, 60 Co spectra were measured in a simple experimen- 
tal setup and compared to simulations. The same setup was thereafter also 
used to study the influence of the source backing by adding copper foils of 
different thickness behind the source. Finally 60 Co (3 spectra were measured 
in strong magnetic fields inside the dilution refrigerator and again compared 
to simulations. All these steps are discussed in detail in the following sections. 

It is important to note that when this work was started the latest available 
version of GEANT4 was GEANT4.7.1. During the writing of this paper a 
new version was released, i.e. GEANT4.8.1. Most of the simulations were then 
redone and compared to the simulations performed with the earlier version. 
Although some small deviations between the two results were found, there 
were no major discrepancies that would change the conclusions of this work. 
Since the comparison of two versions of GEANT4 is not the aim of this paper 
this is not further discussed. For all GEANT4 simulations the low-energy 
electromagnetic processes package was used (j28l ). 



3 GEANT simulations of 60 Co (3 spectrum. 

3.1 Electrons backscattering simulations 

One of the factors that may significantly modify the electron angular distri- 
bution pattern is backscattering of the electrons on the host foil material into 
which the 60 Co ions are diffused, as well as on the layer of solder, on the 
sample holder and on the detector. Therefore, one of the major requirements 
for a good simulation program is a good description of the backscattering of 
electrons. The development of the program was thus started with simulations 
to reproduce data on electron backscattering coefficients already available in 
the literature. Since in the measurements a Si PIN-photodiode was used as a (3 
particle detector, backscattering on silicon was of particular interest. To obtain 
good agreement between our simulations and experimental results available in 
the literature two parameters in GEANT4 were tuned: the CFS parameter 
and the /^-parameter, that were defined above. The electron backscattering 
fractions obtained for different values of these two parameters are listed in 
Table [TJ The simulations were performed for normal incidence of electrons 
with energy E e = 250 keV. The corresponding backscattering fraction calcu- 
lated with the formulas given in Ref. (1291 ) is 12.9 %. Comparing the results 
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of our simulations with this value it may be noted that the value for the f r 
parameter should not exceed 0.02. From the practical point of view, however, 
also the simulation time should be taken into account. Indeed, it was found 
that the simulation time increases with decreasing values for the f r and CFS 
parameters, with the influence of the first one being largest. In view of this 
the f r parameter was then fixed at the value of f r = 0.02. 



In choosing the optimal CFS value for f r = 0.02 one has more freedom. 
Indeed, the calculated backscattered fraction of 12.9 % is determined with an 
absolute precision of only about a few percent (see Ref. (lii!)). Within this 
precision all CFS values in Table [1] give reasonable results for the case with 
f r = 0.02. The CFS value of 10 /xm was then preferred as it also corresponds 
to a reasonable simulation time. This arbitrariness in the choice of the CFS 
parameter value might cause a systematic effect in the determination of the 
(3 asymmetry parameter and the corresponding systematic error on the final 
result of an experiment will have to be determined. 



The optimization of the values for the CFS and f r parameters described here 
was started with GEANT4.7.1. In that version of the toolkit the default value 
of the f r parameter was equal to 0.2. However, as becomes clear from Tabled] 
the adjustment of this parameter is of great importance. In the GEANT4.8.1 
release the default value of f r was set to 0.02 which corresponds to the optimal 
value we also found already in our initial simulations with GEANT4.7.1, as 
well as with the value recommended in Ref. (I24l ). 



In Table [2] our simulation results for the backscattering fraction for two differ- 
ent energies of incident electrons (normal incidence) are compared to results 
obtained by other authors. For these simulations the f r parameter and the 
CFS value were set to 0.02 and 10 /im, respectively. Our results are com- 
pared with the experimental data and simulation results of (jiol ) and with the 
semi-empirical calculations of Ref. (J29i). Good agreement is observed between 
the results of this work and results obtained by these other authors. 



Finally, our simulations showed that the simulated backscattering fractions of 
electrons for large angles of incidence (i.e. above 30° with respect to the normal 
of the detector plane) are larger than the backscattering fractions calculated 
using the equations in Ref. (29) (for the energy dependence) and Ref. (jilj ) (for 
the angular dependence) by a factor of approximately 1.2. 
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Table 1 

Total backscattering fraction (in percent) for different values of the cut-for- 
secondaries CFS and the f r parameter for electrons with kinetic energy E e = 250 
keV at normal incidence. GEANT4.8.1 with the low-energy electromagnetic package 
was used. 







fr 




CFS, /im 


0.2 


0.02 1 


0.002 


10 


8.70(2) 


13.20(3) 


13.20(4) 


5 


9.10(2) 


12.60(3) 


13.50(4) 


1 


10.90(2) 


12.50(3) 


14.30(4) 



1 default GEANT4.8.1 value 

Table 2 

Total backscattering fraction (in percent) of primary electrons with different energies 
E e for normal incidence on to silicon. 



E e , keV 


This work 


Tabata 1 


Berj: 


9 

;er 




(simul.) 




(simul.) 


(expt.) 


250 


13.2 


12.9 


13.2 


14(1) 


500 


11.6 


11.4 


12.3 


14(1) 



1 values calculated with the semi-empirical relation in Ref. II29Q: 

2 fraction of events without total energy absorption (see Ref. d 30f) for details). 

4 Simplified geometry for the 60 Co spectrum measurement and 
simulation 

The next step in the development of the simulation code was to reproduce 
the 60 Co spectrum measured in a simple geometry and with the same detector 
as will be used in real measurements. The setup installed for this purpose is 
shown in Fig. [2j 

4-1 Source preparation and detector configuration 

The source was prepared by drying the 60 Co activity on a 8 /im thick mylar foil 
which was then covered by another mylar foil of the same thickness to prevent 
contamination of surrounding materials. The diameter of the 60 Co spot was 
about 2 mm. The source was attached to a plastic support structure with a 10 
mm diameter hole behind the activity spot to avoid backscattering. This source 
was considered to be a "scattering free" source in all later measurements, but 
the two 8 /im thick mylar foils were of course included in the simulations. A 
500 /im thick Si PIN photodiode with a sensitive area of 81 mm 2 was used as 
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Fig. 2. Schematic view of the test setup for the 60 Co j3 spectrum measurements. 

(3 particle detector. It was mounted on an Al support structure and covered 
by a 1 mm thick copper collimator with 8 mm diameter hole. The detector 
and source were axially aligned with a distance of 50 mm between them. They 
were placed in a vacuum chamber and all measurements were performed at a 
pressure below 1CT 4 torr. 



In all GEANT4 simulations presented in this work the geometry of the setup 
was included in detail. Special care was taken to describe the particle detector 
geometry as precisely as possible. For that purpose HAMAMATSU Photonics 
has provided detailed information on the detector's design (depletion thick- 
ness, front and rear dead layers composition and thickness) which was then 
implemented in the simulation code. 

For comparison with the experimental data the simulated spectra were rescaled 
relative to the experimental ones. The scaling factor was calculated by nor- 
malizing the integral of the simulated spectrum with that of the experimental 
spectrum in the energy region from 150 keV to 300 keV. The same energy 
region was later used in the analysis of the anisotropy of the 60 Co (3 radia- 
tion. The reduced x 2 value for that energy region was calculated taking into 
account the finite statistics of the simulations, which is less than that of the 
experimental data for all simulations shown in this work. 



4-2 Compton background 



The (3 spectrum measured with this simple setup is shown in Fig. [3j The tail 
that extends above the endpoint energy of 60 Co is mainly due to Compton- 
scattered electrons from the 1173 keV and 1332 keV 7 lines that accompany 
the 60 Co (3 decay. The observed (3 spectrum thus contains electrons from both (3 
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Fig. 3. Co (3 spectrum obtained with the "scattering free" source and the Si PIN 
photodiode in the test setup. The regions of the spectrum below and above the 
endpoint energy of 60 Co (i.e. Eq = 318.1 keV) are labeled as 1 and 2 respectively. 
More details are given in the text. 



decays and Compton events. To ensure that the tail is indeed due to Compton 
scattering of the 1173 keV and 1332 keV 7 lines a separate measurement, with 
the detector now covered by a 0.4 mm thick copper plate in which all 60 Co (3 
particles are stopped, has been performed. The result is shown in Fig. HI The 
grey line represents the experimental data while the solid line shows the result 
of the simulations. The detector energy resolution of 7 keV, as determined 
from measurements with a Bi conversion electron source, was included in 
these simulations. 



As can be seen from Fig. HI good agreement with experiment was obtained 
in the energy region of interest (i.e. 150 - 300 keV). Note that at the lowest 
energies the experimental spectrum is disturbed by the tail of the electronic 
noise, which is not included in the simulations. Further, the lowest energy 
part of the spectrum is dominated by (back) scattered and partially absorbed 
(Compton) electrons. It might be that GEANT4 is not perfectly describing 
these processes at such 'low' energies. It is to be noted here that the Compton 
tail is a smooth function of energy over the entire energy region of interest. This 
is very important for the analysis of the anisotropy of the 60 Co (3 radiation 
since the background under the (3 spectrum can cause some uncertainty in 
the determination of the absolute value of the (3 anisotropy and thus in the 
determination of the (3 asymmetry parameter A. 
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E, keV 

Fig. 4. Comparison of the experimental and simulated Compton background from 
60 Co. Grey dotted line: experimental data (the width of the line is due to the la 
statistical error bar); black solid line: GEANT4.7.1 simulations. The reduced x 2 f° r 
the simulations when compared to the data in the energy range 150 - 300 keV is 
also indicated. 

4-3 (3 spectrum 

The next step after simulating the Compton background was to reproduce 
the full j3 spectrum of 60 Co shown in Fig. [3j The result is presented in Fig. [5j 
Possible reasons for the discrepancy in the region above 300 keV are discussed 
below. 

4-4 Effect of backs cattering 

To estimate the influence of the source backing a series of measurements was 
performed. In these the 60 Co j3 spectrum was measured with copper foils of 
different thicknesses (i.e. 10, 20, 50, 100 and 400 /mi) installed behind the 
scattering-free source. Simulations for these cases were then performed as well 
and compared to the experimental data. The results are shown in Fig. [6l The 
experimental data (panel " a" ) are normalized to the counting time. The sim- 
ulated spectra shown in panel "b" are again rescaled to the corresponding 
experimental spectra by normalizing the integrals of the simulated and exper- 
imental spectra in the energy region from 150 keV to 300 keV, as mentioned 
above. 

From the data shown in Fig. [6] the general effect of increasing the thickness of 
the backing is clearly visible: thicker backings cause a larger intensity in the 
low-energy part of the /3 spectrum (caused by an increase in the backscatter- 
ing from the source backing) as well as a larger contribution from Compton 
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Fig. 5. (a) Comparison of the experimental spectrum of 60 Co of Fig. [3] with the 
simulated one. Grey dotted line: experimental data (error bars are smaller than the 
size of the data points). Black solid line: GEANT4.8.1 simulations with f r = 0.02 
and CFS = 10 //m. (b) Relative difference between simulated and experimental 
spectra. In the region of interest the difference is within 2 to 3 %. It is larger near 
the (3 spectrum endpoint where the relative contribution of Compton background 
is more important. The value x 2 l v = 1-3 for 67 degrees of freedom (which is within 
the 95% confidence level) refers to the energy region of interest, i.e. from 150 keV 
to 300 keV. 



events. The backscattered fraction between 50 and 150 keV is clearly satu- 
rating already between 50 and 400 /im, a feature that is reproduced by the 
simulations (see Figs. [6]^" a" ) and E]("b")). A direct comparison of experimen- 
tal data and simulations for Cu backings of 10 fim and 400 /im thickness is 
shown in Figs. [6]("c") and [B("d"). It is seen that the simulations do very well 
between 50 and 300 keV for all backing thicknesses, in particular in the region 
between 50 and 150 keV which is directly sensitive to backscattering. 

Above 300 keV (i.e. in the region that is dominated by Compton electrons) the 
data/simulation ratio is found to be about 1.5 (see e.g. Fig. [5^ and Figs. 
and [6]i), independent of backing thickness. Simulations showed that for the 
setup in Fig. [2] that was used to measure the spectra shown in Figs. [5] and EJ 
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the source foil and the detector contributed about equally to the number of 
detected Compton electrons. The difference of about 50% between the number 
of Compton electrons in the data and the simulations then suggests that either 
the amount of Compton scattering of 7 rays, or the description of multiple 
scattering of the Compton electrons, before they reach the detector, is not 
fully adequate (see also Ref. 

However, this discrepancy between the data and simulations with respect to 
the number of Compton electrons will not have an important effect on the 
analysis of the 60 Co (3 anisotropy if the Compton background can be properly 
subtracted. In order to do so one has to know (1): the exact shape of the 
Compton background in the energy region of the 60 Co (3 spectrum and (2): 
the amplitude of this Compton background. As was shown above (see Fig. 
Uj) GEANT4 reliably describes the Compton background down to energies of 
about 100 keV, thus providing the shape of it. Further, knowing the shape of 
the Compton background its amplitude can be obtained from the part of the 
experimental spectrum above the /3 endpoint. Applying this procedure to the 
data shown in Fig. causes the reduced \ 2 t° increase from 1.3 to 1.9. The 
subtraction is thus not perfect, and for analysis of real data an error related 
to this subtraction procedure will be included. 



5 Simulation of the 60 Co spectrum in the real setup 

As a final step in the development of the simulation code that should take 
into account the effects that modify the experimental (3 spectra, the 60 Co 
13 spectrum was measured in the actual setup, this time also with non-zero 
magnetic field, and compared to simulations for these new and more delicate 
conditions. Three additional factors now had to be considered: 

• a realistic description of the 60 Co source; 

• the description of the propagation of the f3 particles in the magnetic field; 

• the anisotropy of the (3 radiation emitted from a polarized sample. 

These issues are considered in detail below. 



5. 1 Description of the source 

The 60 Co source was prepared by diffusing the 60 Co activity into a 20 /zm thick 
copper foil that was maintained at a temperature of 850°C for 5 minutes. This 
procedure allowed to prepare a source with the radioactive 60 Co nuclei being 
contained in an about 3 /zm thin layer below the surface of the copper foil. 
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Fig. 6. Experimental (panel A) and simulated (panel B) 60 Co [5 spectra for the 
"scattering-free" source (see text) backed by copper foils with different thicknesses. 
The two bottom panels compare experimental data and simulations for a Cu backing 
of 10 //m (panel C) and 400 //m (panel D). Note that for the spectrum obtained 
with a 400 j^m backing, where the scattered fraction is largest, the reduced x 2 is 
slightly worse. 



The diffusion of cobalt in copper was studied in detail in Ref. (]32l ). Using the 
diffusion coefficient listed there for the temperature of 850 °C at which our 
source was produced, we calculated the diffusion depth profile of the cobalt 
ions for this temperature and the diffusion time that was used. This profile 
was then implemented in the simulations. 



5.2 Description of the propagation of the (3 particles in the magnetic field 



For the description of the magnetic field the field map in radial and axial 
directions with a step of 1 mm, provided by the manufacturer of the mag- 
net (Oxford Instruments Ltd.), was implemented in the simulation code. In 
the GEANT4 toolkit Runge-Kutta integration is used to compute the motion 
of a charged particle in a general field (1331 ) and different methods of inte- 
gration are provided. The specific method (or "stepper") can be chosen by 
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the user, with the preferred choice being determined by the particular field 
conditions. In this work four different steppers provided by GEANT4 were 
compared: GJ^HelixExplicitEuler, G4HelixImplicitEuler, G4HelixSimpleRunge 
and G4ClassicalRK4 (default stepper). The first three are specifically de- 
signed for the description of particle propagation in a purely magnetic field 
while the last one is a default GEANT4 stepper that can be used for general 
electromagnetic fields. The GEANT4.8.1 simulated spectra of 60 Co obtained 
with the four different steppers were found to be very similar to each other. 
This is at variance with our earlier simulations performed with GEANT4.7.1 
where the spectra simulated with different steppers were clearly differing. The 
best agreement with the experimental spectrum was then obtained for the 
G4HelixExplicitEuler. This same stepper was also used in all later simulations 
with GEANT4.8.1. 



Typical trajectories of electrons with energies of 100 keV and 300 keV in a 13 
T magnetic field are shown in Fig [3 The positions of the decaying 60 Co nuclei 
emitting the electrons are chosen randomly in a circle with a diameter of 1 
mm. The thickness of the copper host foil is 20 fim and the depth at which the 
decaying nuclei are sitting is randomized according to the depth distribution 
profile discussed above. 



The different helix radii of electrons spiraling in the magnetic field are caused 
by the different energies and angles at which electrons leave the foil. Because 
of the very large magnetic field strength the electrons are strongly focused, 
leading to an effective detection solid angle of almost 2n. Further, due to the 
divergence of the magnetic field lines towards the detector and the principle 
of adiabatic invariance of the magnetic flux (1341 ). i.e. p 2 ± / B=const. (with p± 
the momentum component perpendicular to the magnetic field B), most of 
the transverse momentum of the electrons is transferred into longitudinal mo- 
mentum. Therefore, no electrons arrive at the detector at angles of incidence 
with respect to the normal on the detector surface larger than about 20°. 



An interesting feature may finally be seen in the bottom panel of Fig. i.e. 
the backscattering of an electron on the sample holder (trajectory extending 
in the negative Z direction). The sample holder has a 9 mm deep hole with 
a diameter of 6 mm and a rough surface (screw thread) behind the foil to 
avoid such backscattered events, but these can of course still occur from time 
to time on the bottom surface of this hole (see Fig. [1]). Simulations showed, 
however, that this happens for only about one permille of the electrons that 
finally arrive in the detector. 
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Fig. 7. Trajectories of electrons with energies of 100 keV (upper panel) and 300 keV 
(bottom panel) in a magnetic field of 13 T. X-Z plane; the magnetic field is along 
the Z-axis. 

5.3 Anisotropy of the (3 radiation 



To fully implement the GEANT4 based method described above for the anal- 
ysis of experimental data, also the anisotropy of the (3 radiation finally had to 
be included in the simulations. The anisotropy of the electrons emitted from 
60 Co was calculated using Eq. [1] assuming the absence of any kind of exotic 
physics (i.e. assuming a pure V — A interaction with maximal parity viola- 
tion and no time-reversal violation). All terms in Eq. [T] are then known and 
the angular distribution of the 60 Co (3 decay electrons as a function of their 
emission angle 9 can be implemented in the simulation code. 

In Fig. M the simulated isotropic (i.e. unpolarized nuclei, panel "a") and 
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anisotropic (polarized nuclei, panel "b") spectra are compared with experi- 
mental data obtained with 60 Co. For the experimental spectra a pile-up prob- 
ability of 0.18 % was obtained from the shaping time and count rate. This 
value was then used in the simulations. In the analysis of future measure- 
ments we will assign a generous error (i.e. of 50 %) to this estimated pile-up 
rate. Further, as we are always measuring at not too high counting rates (typi- 
cally about 1 kHz or lower) the pile-up correction will not induce a large effect 
in the simulated spectra and will therefore not have a significant effect on the 
final result. 



In Ref. ( 124j ) it was shown that simulation results do not change anymore when 
values for the cut-for-secondaries parameter CFS smaller than 1 fim are used. 
For such small values simulation times become exceedingly large, however. 
Simulations were therefore performed for CFS = 1 /im and 10 /im. As no 
difference could be observed within the statistical errors of these simulations 
the value of the CFS parameter was set to 10 /im so as to obtain reasonable 
simulation times. 



The experimental spectra in Fig. [H] are "distorted" in the energy region below 
150 keV by conversion electrons from the 57 Co activity that was diffused in 
the same foil as 60 Co, for thermometry. Good agreement between experimental 
data and simulations is again obtained in the energy region from 150 to 300 
keV, in spite of the fact that the amplitude of the Compton tail above the spec- 
trum endpoint is again too low in the simulated spectra. In Fig.Oit was shown 
already that the Compton fraction depends on the source conditions (as e.g. 
the thickness of the backing), but there may be more parameters determining 
this. Simulations showed that for the setup with the 13 T magnetic field that 
was used here about 99 % of the detected Compton electrons originate from 
the Cu source foil (compared to an about equal number of Compton electrons 
being produced in the source foil and in the detector when no magnetic field 
is present; cf. end of section 4). This can be easily understood. A Compton 
electron created in the source will be focussed directly on to the detector by 
the strong magnetic field, similar to the (3 decay electrons. However, in order 
to have Compton scattering occur in the detector, 7 rays from the sample first 
have to reach the detector, which has a very small (geometrical) solid angle 
(order 10~ 4 ). The remaining 1 % of the detected Compton electrons mainly 
originate from the Cu sample holder. Compton electrons created elsewhere in 
the setup were found to have a negligibly small probability to make it to the 
detector. 

Simulations also showed that the presence of the weak 54 Mn thermometer 
(with a 835 keV 7 ray) caused a negligible amount of Compton electrons 
reaching the detector. This can be easily understood from the geometry of 
the setup as shown in Fig. [TJ most Compton electrons from gamma rays of 
the 54 Mn source are indeed absorbed in the part of the sample holder that is 
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situated between the 54 Mn source and the Si PIN diode detector. 

The 57 Co conversion electron peaks could be well reproduced in the simula- 
tions. But to make these cause a significant increase in the amount of counts 
above the spectrum endpoint in the simulations, a much too large pile-up per- 
centage had to be assumed (i.e. of the order of 1 % or more, compared to an 
expected pile-up rate of about 0.2 % given the count rates and amplifier shap- 
ing times). One is therefore again lead to conclude that the counts above the 
beta spectrum endpoint most probably are Compton electrons related to the 
gamma rays of 60 Co. Then, either the amount of Compton scattering of 7 rays 
in the backing, or the description of multiple scattering so that they can reach 
the detector, is not fully adequate. It was found that when assuming the ac- 
tivity to be distributed throughout the entire foil (which would have required 
a much longer diffusion time than was actually used), leading to an increased 
amount of Compton scattering in the backing, the amount of Compton counts 
in the spectrum relative to the number of detected (3 particles increases only 
by about 20 %. This is less than half of the 50 % increase that is needed. This 
could indicate that it is the description of multiple scattering of the Compton 
electrons, so that they can reach the detector, that is not fully adequate. It 
is not clear yet why the simulation code does not describe this correctly and 
of course this will have to be taken into account as a systematic error in the 
analysis of experimental data with this code. 



6 Conclusion 



A GEANT4 based Monte-Carlo routine was developed and optimized to sim- 
ulate the spectrum shape for (3 particles emitted by unpolarized as well as 
polarized nuclei in LTNO experiments and to extract the f3 asymmetry pa- 
rameter using Eq. [31 The code accounts for the different effects that modify 
the emitted (3 spectrum, i.e. (back) scattering, magnetic field effects, etc. 
Good agreement between simulations and experimental data in the energy re- 
gion of interest is found for a 500 /im Si PIN diode detector with the relative 
difference always b eing less than about 2%. Similar results were obtained by 



other authors (j22|)-(|25|). 



Due to the modularity of this simulation code and the object-oriented struc- 
ture of the GEANT4 toolkit this method can be easily adopted for similar 
types of measurements with a different geometry and a different isotope stud- 
ied. In particular it is now being implemented for a series of precision mea- 
surements of the (3 asymmetry parameter of selected isotopes that is being 
performed both in Leuven and at the NICOLE low temperature nuclear ori- 



entation setup at ISOLDE/CERN (jilt: Il6l). The code can in principle rather 



easily be generalized to other be used also for experimental setups using other 
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polarization techniques. 
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Fig. 8. Comparison of the experimental and simulated Co spectrum for the actual 
experimental setup and a magnetic field B = 13 T. Panel a: "warm" data (isotropic 
emission, unpolarized nuclei; x 2 l v = 1-2 for 40 degrees of freedom in the region of 
interest from 150 to 300 keV); panel b: "cold" data (anisotropic emission, polarized 
nuclei; x 2 l v = 0-8 for 40 degrees of freedom in the region from 150 to 300 keV). The 
peaks between 100 keV and 150 keV in the experimental (5 spectrum are conversion 
electrons from 57 Co. The Compton tail in the experimental spectra also includes a 
small amount of Compton background from the 54 Mn thermometer source that was 
also present on the sample holder. Note that the 57 Co and 54 Mn sources were not 
included in the simulations. 
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